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Abstract 

Incremental  collapse  analysis  invariably  involves  simulation  of  failure  mecha¬ 
nisms  of  'distributed*  or  highly  ' localized *  nature. 

The  report  describes  current  work  on  the  progressive  failure  behavior  of  con-  _ 
crete.  After  a  brief  review  of  experimental  failure  observations  in  tension  and  >n 
compression  the  issues  of  ' normality *  and  ' characteristic  length*  are  discussed 
when  describing  localized  changes  of  the  internal  structure  during  progressive  fail¬ 
ure.  To  this  end  the  internal  microstructure  of  concrete  is  considered  to  form  a 
two-phase  particulate  composite  made  of  elastic  aggregate  embedded  in  inelastic 
mortar.  The  paper  examines  the  degradation  of  macroscopic  stiffness  and  strength 
with  the  aid  of  numerical  simulation  of  the  mortar-aggregate  composite.  The  com¬ 
putational  results  provide  new  insight  into  the  distributed  and  localized  nature  of 
tensile  and  compressive  failure  processes  for  comparison  with  experimental  obser¬ 
vations  on  concrete  composites. 
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In  view  of  the  importance  of  reliable  safety  assessments  there  is  an  increasing  demand 
for  accurate  failure  predictions  which  trace  the  response  behavior  beyond  the  frame¬ 
work  of  well-posed  initial  boundary  value  problems.  Progressive  de-stabilization  of 
the  spatial  and  temporal  discretization  invariably  leads  to  loss  of  single-valued  solu¬ 
tion  regimes.  It  remains  to  construct  computational  strategies  which  not  only  survive 
but  also  reproduce  the  actual  failure  progression. 

The  computational  difficulties  are  compounded  by  the  fact  that  peak  and  post¬ 
peak  response  phenomena  are  really  the  manifestations  of  dramatic  changes  of  the 
microstructure  which  can  no  longer  be  considered  material  properties.  One  of  the 
main  difficulties  is  therefore  the  extraction  of  fundamental  fracture  properties  from 
servo-controlled  laboratory  experiments  on  specimens  which  are  tested  beyond  peak 
into  the  post-peak  regime.  Backward  identification  of  the  governing  fracture  param¬ 
eters,  such  as  the  characteristic  width  of  the  localized  failure  band  and  the  dilatancy 
during  progressive  failure,  involves  the  analysis  of  discontinuities  and  is  one  of  the 
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pressing  issues  of  failure  mechanics,  see  recent  review  by  Bazant  [l]. 

In  the  case  of  concrete  and  other  cementitious  composites  very  different  failure 
modes  can  be  mobilized.  They  depend  primarily  on  the  type  of  loading  and  range 
from  very  localized  mechanisms  in  tension  to  very  distributed  failure  in  confined  com¬ 
pression.  The  success  of  computational  failure  simulations  depends  critically  on  the 
macroscopic  constitutive  model.  Aside  from  proper  initiation  conditions  realistic  fail¬ 
ure  formulations  have  to  capture  fracture  propagation  and  thus  the  drastic  transition 
from  very  brittle  to  very  ductile  fracture  as  a  fuction  of  confinement,  see  Hegemier 
and  Read  [2]  as  well  as  Horii  and  Nemat-Nasser  [3j. 

In  the  first  part,  the  paper  summarizes  experimental  observations  of  triaxial  failure 
in  tension  and  compression.  In  the  second  part,  current  constitutive  developments  are 
briefly  discussed  with  respect  to  their  ability  to  predict  failure.  In  the  third  and  last 
part,  the  question  is  examined,  whether  the  heterogeneous  nature  of  concrete  does 
affect  initiation  as  well  as  propagation  of  failure  in  tension  and  compression.  Thereby 
it  is  understood,  that  the  size  of  the  largest  aggregate  introduces  an  internal  mate¬ 
rial  scale  which  controls  the  distribution  of  initial  microdefects  and  thus  the  location 
of  initial  failure  mechanisms.  The  question  arises,  to  what  evtent  does  the  internal 
micro-structure  influence  the  critical  fracture  path  in  tension  and  compression  dur¬ 
ing  coalescence  of  microdefects  into  the  critical  macrodefect?  In  broader  terms,  we 
ask  whether  the  macroscopic  failure  mode  depends  on  the  internal  micro-stucture  and 
the  “critical  failure  path "  in  the  particulate  composite  and  whether  the  macroscopic 
degradation  of  strength  and  the  dilatancy  are  affected  altogether? 

1.1  Experimental  Observations  of  Concrete  Failure 

The  failure  behavior  of  concrete  is  governed  by  complex  degradation  processes  within 
the  matrix-aggregate  composite. 

In  direct  tension,  we  observe  a  very  brittle  separation  process  which  results  in 
the  formation  of  highly  localized  tensile  cracks.  The  direct  tension  experiments  by 
Reinhardt  [4],  Gopalaratnam  and  Shah  [5]  and  Hurlbut  [6]  on  cylindrical  specimens  of 
different  height  showed  that  the  formation  of  macroscopic  discontinuities  is  attributed 
to  a  surface-dominated  failure  process  in  the  spirit  of  fracture  mechanics.  In  fact,  the 
results  of  the  direct  tension  experiments  depicted  in  Fig.  1  illustrate  the  brittle  nature 
of  concrete  with  a  well  defined  and  stable  post-peak  regime.  They  exhibit  a  surpris¬ 
ing  amount  of  residual  tensile  strength  due  to  ‘teething’  with  relatively  little  stiffness 
degradation  before  and  after  peak.  In  sum,  tensile  cracking  is  a  fracture  phenomenon  in 
spite  of  the  tortuous  character  of  the  uack  surfaces  in  the  matrix-aggregate  composite. 

At  low  confined  compression,  we  observe  a  more  gradual  degradation  process  of 
strength  which  is  accompanied  by  highly  dilatant  behavior  in  the  lateral  direction. 
The  uniaxial  compression  experiments  by  von  Mier  [7]  on  prismatic  specimens  of  dif¬ 
ferent  height  suggest  that  compressive  failure  is  also  dominated  by  fracture  processes 
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which  are  however  more  distributed  as  we  increase  the  level  of  lateral  confinement. 
F ig. 2  illustrates  the  large  differences  in  the  load-deformation  behavior  of  cylindrical 
specimens  subject  to  different  levels  of  lateral  confinement.  The  experimental  results, 
which  are  extracted  from  the  comprehensive  test  series  by  Smith  [9],  exhibit  the  rapid 
increase  of  strength  and  ductility  with  virtually  intact  elastic  unloading  behavior.  In 
fact,  there  is  a  transition  point  TP  of  brittle-ductile  fracture  which  separates  softening 
of  the  macroscopic  strength  from  highly  ductile  failure  regimes  in  which  virtually  no 
degradation  of  strength  is  observed. 

The  schematic  plot  in  Fig.3  delineates  the  different  failure  modes  in  the  P.cr.dulic 
sUcaS  plane.  The  narrow  softening  zone  extends  from  tension  into  the  low  confined 
compression  regime  up  to  the  transition  point  of  brittle-ductile  fracture.  The  failure 
regions  beyond  the  transition  point  are  very  ductile  and  engulf  a  wide  region  of  hard¬ 
ening  behavior  in  triaxial  compression  which  are  compactive  near  the  hydrostat  and 
which  turn  gradually  dilatant  in  the  vicinity  of  the  failure  envelope.  These  observations 
are  based  on  pertinent  triaxial  load  history  experiments  which  have  been  computerized 
by  Stankowski  [8]  in  the  “ triaxial  data  base  manual *  for  concrete  comprising  stress- 
strain  histories  of  136  experiments.  The  global  failure  mechanisms  below  the  transition 
point  differ  strongly  when  the  response  in  triaxial  compression  is  compared  with  that 
in  triaxial  extension.  In  fact,  cylindrical  fracture  in  the  form  of  * orthogonal ’  cleavage 
is  responsible  for  rather  ductile  softening  on  the  macroscopic  level.  This  failure  in  tri¬ 
axial  compression  is  primarily  due  to  the  excessive  lateral  tensile  strains  in  the  radial 
and  circumferential  directions  when  confined  cylindrical  test  specimens  are  subjected 
to  increasing  axial  compression.  On  the  other  hand,  “ planar  "  separation  is  responsible 
for  very  brittle  softening  on  the  macroscopic  level  due  to  excessive  axial  tensile  strains 
when  cylindrical  specimens  are  subject  to  increasing  lateral  confinement  in  triaxial 
extension. 

This  overview  of  triaxial  deformation  and  failure  mechai  isms  sets  the  stage  for 
physically  motivated  concrete  model  which  should  reproduce  the  following  features  in 
the  pre-  and  post-peak  regimes: 

•  There  exists  a  transition  point  of  brittle-ductile  fracture  which  separates  soften¬ 
ing  from  hardening  behavior, 

•  The  degradation  of  strength  in  tension  as  well  as  in  compression  is  caused  by 
fracture  processes  of  the  internal  composite  structure  resulting  in  macroscopic 
softening, 

•  There  are  two  distinct  modes  of  cylindrical  and  planar  failure  which  are  asso¬ 
ciated  with  orthogonal  and  planar  cleavage  in  triaxial  compression  and  triaxial 
extension, 

•  The  lateral  deformation  behavior  is  characterized  by  loss  of  normality.  Thus 
non-associated  flow  determines  the  amount  of  inelastic  dilatancy  and  also  the 
triaxial  strength  when  strain  control  prevails. 
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In  sum,  the  pressure-sensitive  failure  envelope  differs  markedly  between  loading 
paths  in  triaxial  compression  and  triaxiai  extension.  Moreover,  hardening  as  well  as 
softening  response  functions  depend  strongly  on  the  lateral  confinement.  Thus  the 
traditional  plasticity  concept  of  a  ‘universal*  equivalent  stress-strain  relation  has  to  be 
extended  in  order  to  account  for  the  large  increase  of  ductility  with  increasing  con¬ 
finement.  Finally,  the  nonlinear  compactive  behavior  in  the  vicinity  of  the  hydrostat 
supports  the  notion  of  a  cap.  The  inelastic  dilatancy  in  the  dilative  cone  and  the 
compactive  cap  regions  however  requires  direct  control  in  the  form  of  non-associated 
flow  separately  from  the  yield  condition. 

1.2  Normality  of  Inelastic  Dilatancy 

The  non-associated  formulation  determines  primarily  the  amount  of  inelastic  dilatancy 
and  thus  the  lateral  deformation  behavior  in  stress-controlled  environments.  On  the 
other  hand,  in  strain-controlled  environments  it  is  important  to  realize  that  normality 
severely  overpredicts  the  triaxial  strength  because  of  the  excessive  dilatation  which  is 
converted  into  excessive  lateral  confinement.  Therefore,  inelastic  dilatancy  control  is 
not  merely  benefiting  the  lateral  deformation  characteristics,  it  determines  the  strength 
and  thus  the  ductility  of  the  response  predictions  in  strain-controlled  environments. 

Fig. 4  taken  from  the  recent  experiments  by  S.  Smith  [9]  illustrates  the  directions  of 
the  inelastic  strain  rates  which  are  superimposed  on  the  corresponding  stress  path  of 
the  triaxial  compression  tests  at  different  confinement.  Assuming  negligible  degrada¬ 
tion  of  the  elastic  properties  the  total  deformation  measurements  can  be  decomposed 
into  the  sum  of  elastic  and  inelastic  components.  The  angle  of  the  inelastic  strain  rate 
component  with  the  hydrostat  provides  a  measure  of  dilatancy  which  is  invariably 
smaller  than  the  corresponding  friction  angle,  i.e.  the  normal  to  the  current  loading 
surface.  This  experimental  observations  agrees  with  the  general  opinion  that  frictional 
materials  are  non-associative  and  involve  therefore  non-symmetric  constitutive  opera¬ 
tors.  Another  feature  of  this  loss  of  normality  is  depicted  in  Figs. 5  and  6  where  the 
results  of  the  “wedge  test’  are  shown.  This  experiment  was  carried  out  in  response  to 
the  recent  discussion  by  Sandler  and  Rubin  [11]  who  demonstrated  analytically  that 
loading  into  the  wedge  region  is  synonymous  to  negative  energy  dissipation  and  thus 
leads  to  instability.  Fig. 6  verifies  this  argument  with  respect  to  the  loss  of  positive 
plastic  work  increment,  i.e  d7Wp  =  /  aedv  <  0.  However,  this  is  only  a  sufficient  con¬ 
dition  and  not  a  necessary  one.  In  fact,  the  specimen  remained  fully  intact  under  load 
control  inspite  of  the  negative  incremental  plastic  work. 

In  view  of  the  prevalent  displacement  setting  of  the  majority  of  finite  element  pro¬ 
grams  normality  or  rather  the  loss  of  normality,  is  an  extremely  important  feature 
which  can  alter  totally  the  triaxial  response  predictions.  The  dramatic  effect  of  asso¬ 
ciated  versus  non-associated  flow  is  illustrated  with  a  model  problem  taken  from  the 
recent  thesis  by  Pramono  [12].  In  that  case  the  Leon  failure  model  [13]  is  subjected  to 
a  strain  path  extracted  from  a  triaxial  compression  experiment  at  aT  =  200  psi  con¬ 
finement.  Assuming  in  both  cases  perfectly  plastic  behavior  the  response  predictions 
of  associated  versus  non-associated  flow  are  totally  different,  when  the  plastic  poten- 
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tial  of  the  non  -associated  formulation  is  governed  by  the  pressure  insensitive  Tresca 
criterion.  The  response  predictions  are  shown  in  Ftg.7  for  the  two  models.  After  the 
failure  envelope  has  been  reached,  the  stress  path  of  the  associated  flow  rule  moves 
up  the  failure  envelope  until  a  point  is  reached,  where  the  prescribed  strain  increment 
is  equal  and  coaxial  with  the  plastic  strain  increment.  On  the  stress-strain  level  this 
results  in  an  excessive  axial  strengt1.  prediction  with  unrealistically  ductile  behavior 
in  the  form  of  perfectly  plastic  flow.  Conversely,  the  non-associated  flow  rule  moves 
the  stress  path  down  the  failure  envelope  into  the  tension  region.  On  the  stress-strain 
level  this  results  in  a  reasonable  prediction  of  axial  strength  with  post-peak  behavior 
which  is  too  brittle  and  exhibits  too  rapid  degradation  of  the  residual  strength. 

These  totally  opposite  response  predictions  demonstrate  the  fundamental  role  of  the 
dilatancy  on  the  overall  strength  and  ductility  behavior  of  pressure  sensitive  models. 


2  Current  Constitutive  Issues 

Current  concrete  models  are  either  based  on  nonlinear  elasticity,  or  on  alternative 
plasticity,  continuous  damage  and  fracture  mechanics  theories  respectively.  The  large 
number  of  recent  publications  are  thus  of  little  surprise,  ([14],  [15],  [16],  [17],  [18],  [19], 
[20],  [21],  [22],  [23],  [24],  [25])  since  none  of  these  frameworks  adequately  descibes  the 
entire  spectrum  of  deformation  and  failure  behavior  of  concrete.  In  view  of  the  severe 
shortcommings  of  most  formulations,  it  is  imparative  to  verify  these  phenomenological 
constitutive  models  by  micromechanical  investigations,  partly  because  of  the  limited 
number  of  experiments  which  are  available  for  calibration  and  verification  [8]. 

In  the  context  of  micromechanics  it  is  instructive  to  recall  some  major  results  of 
composite  mechanics,  see  ref.  by  Hashin  [26].  In  the  case  of  “ isotropic  composites’ 
made  of  a  matrix  and  much  stiffer  and  stronger  spherical  particles  the  effective  com¬ 
posite  properties  depend  primarily  on  the  volume  fractions  of  the  constituents.  While 
the  elastic  stiffness  properties  can  be  increased  up  to  four  or  five  times  the  matrix 
stiffness,  the  strength  of  such  composites  is  only  of  the  order  of  the  matrix  strength 
and  may  be  even  lower.  The  stiffness  and  strength  of  the  composite  increases  when 
elongated  particles  are  used,  however,  at  the  price  of  lower  attainable  volume  fractions. 
For  particulate  composites  the  fundamental  postulate  of  heterogeneous  media  states: 
“ The  stress  and  strain  fields  in  a  large  heterogeneous  body  subjected  to  homogeneous 
boundary  conditions  are  statistically  homogeneous,  except  in  a  boundary  layer  near  the 
external  surface ”.  The  formulation  of  effective  elastic  properties  of  composite  materi¬ 
als  centers  around  the  concepts  of  “ Representative  Volume  Elements’  and  “ Statistical 
Homogeneity’.  The  notion  of  an  equivalent  continuum  involves  three  characteristic 
length  scales, 

•  the  micro  -scale,  which  is  of  the  size  of  characteristic  particle  diameters, 

•  the  meso(mini)-sca!e,  which  corresponds  to  the  size  of  homogeneized  representa¬ 
tive  volume  elements,  and 

•  the  macro-scale,  which  is  of  the  dimension  of  the  composite  material  body. 
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In  the  following  we  examine  tensile  and  compressive  failure  mechanisms  on  the 
meso-level.  The  geometry  of  the  two-phase  composite  body  is  in  this  case  assumed  to 
be  of  the  same  size  as  the  representative  volume  element.  For  concrete  this  dimension  is 
usually  at  least  six  times  the  maximum  aggregate  size  for  which  statistical  homogeneity 
may  be  adopted.  The  main  intent  is  to  gain  new  insight  into  the  evolution  of  the 
underlying  failure  processes  in  particle  composites  for  which  quantitative  results  are 
extremely  difficult  to  obtain  from  macroscopic  experiments.  The  study  provides  aiso  a 
testbed  for  probing  different  constitutive  hypotheses  such  as  Mode  I  controlled  failure 
mechanisms  and  their  effects  on  the  macroscopic  behavior. 


3  Numerical  Simulation  of  Particle  Composite 

The  numerical  simulation  of  the  two-phase  mortar-aggregate  composite  is  approached 
on  two  levels:  First,  the  geometry  and  topology  of  the  phase  interface  geometry  is 
generated  in  random  manner  in  order  to  define  the  configuration  of  representative 
composite  bodies.  Second,  the  response  behavior  is  evaluated  deterministically  when 
the  composite  samples  are  subjected  to  prescribed  load  histories  probing  elementary 
constitutive  assumptions  for  the  physical  properties  of  the  iwo  constituents. 

3.1  Geometry  and  Topology  of  Mortar- Aggregate  Phases 

Concrete  is  a  multiphase  material  comprised  of  well-graded  aggregate,  hydrated  ce¬ 
ment  paste  and  voids.  Since  not  every  particle  or  void  can  be  considered  explicitely, 
the  numerical  simulation  is  confined  to  a  two-phase  matrix-particle  composite  in  which 
only  the  larger  inclusions  are  accounted  for.  Opposite  to  the  deterministic  genera¬ 
tion  of  Numerical  Concrete ’  based  on  morphological  descriptions  by  Roelfstra  and 
Wittmann  [27]  the  mortar-aggregate  composite  is  here  described  in  terms  of  polygo¬ 
nal  approximations  of  crushed  angular  aggregate  particles  which  are  embedded  in  the 
mortar  matrix  in  a  random  manner.  The  polygonal  description  was  adopted  earlier  by 
Cundall  et  al.  [28]  for  the  idealization  of  granular  solids  with  rigid  particles  and  sliding 
interfaces,  also  known  under  the  name  of  the  “ Distinct  Element  Method ".  In  our  case, 
the  composite  body  is  subdivided  into  Voronoi-Polygons  [29]  centered  around  poles 
which  are  randomly  distributed  trough  Monte  Carlo  perturbation  of  a  regular  grid. 

The  interface  geometry  of  the  mortar-aggregate  composite  is  generated  as  follows: 
A  systematic  decision  tree  identifies  “non-contiguous’  aggregate  polygons  by  shrinking 
the  size  of  the  original  polygons  in  a  random,  non-affine  manner.  Fig. 8  illustrates  the 
composite  geometry  generated  by  this  non-contiguous  approach  allowing  direct  control 
of  the  aggregate  area  fraction.  In  the  subsequent  simulation  an  aggregate  area  fraction 
of  50%  is  used  throughout. 

The  gradation  of  the  composite  is  characterized  by  the  cumulative  size  distribution 
of  aggregate  particles.  The  sieveline  in  Fig. 9  illustrates  the  aggregate  distribution 
generated  by  the  non-contiguous  inclusion  strategy,  whereby  smoother  distributions 
could  have  been  attained  by  refining  the  grid  and  appropriate  spacing.  The  resulting 
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size  distribution  is  certainly  satisfactory  in  view  of  the  limitations  inherent  in  the  two- 
dimensional  idealization  of  the  three-dimensional  microstructure. 

The  two-dimensional  finite  element  discretization  with  2352  DOF  is  shown  in 
.Ft'y.lO  which  also  depicts  the  idealizations  of  the  two  constituents.  Each  polygon 
is  subdivided  into  triangular  element  patches  made  of  four  or  nine  constant  strain 
triangles  which  span  between  two  adjacent  vertices  and  the  centroid  of  the  polygon. 
The  complex  idealization  illustrates  the  need  for  automatic  pre-  and  post-processing 
with  renumbering  and  band  optimization.  We  also  note  the  rapid  growth  of  elements 
and  number  of  degrees  of  freedom,  DOF,  when  the  polygonization  is  refined  for  better 
gradation.  In  fact,  the  more  appropriate  three-dimensional  simulation  of  the  compos¬ 
ite  is  simply  beyond  current  computational  means,  although  conceptually,  it  does  not 
require  development  of  new  features. 

3.2  Material  Description  of  Constituents 

The  main  objective  of  this  study  is  to  evaluate  the  influence  of  the  composite  action  on 
the  failure  mechanisms  in  concrete.  Therefore,  only  elementary  constitutive  assump¬ 
tions  are  adopted  for  each  of  the  constituents  in  order  to  quantify  the  composite  effect 
on  the  macroscopic  manifestations  of  localized  failure. 

In  the  computational  failure  study  the  aggregate-phase  is  assumed  to  remain  in 
the  linear  elastic  response  regime.  This  constitutive  postulate  is  fairly  realistic  for  low 
and  medium  strength  concretes  with  stiffness  and  strength  properties  of  the  aggregate 
particles  which  are  substantially  greater  than  those  of  the  mortar. 

The  scope  of  the  current  composite  study  is  limited  to  an  elastic-plastic  Drucker - 
Prager  formulation  of  the  mortar  matrix  with  decohesive  softening  in  the  post-peak 
regime.  “ Normality *  is  assumed  in  order  to  quantify  the  loss  of  associated  flow  on  the 
macroscopic  level.  Thus,  the  friction  angle  is  kept  constant  during  plastic  flow  and 
coincides  with  a  dilatancy  angle  of  $  =  30®  for  the  matrix  constituent.  Furthermore, 
no  characteristic  length  is  introduced  to  regularize  the  softening  computations  and  to 
minimize  mesh  size  dependence  [30],  [31].  In  fact,  the  composite  argument  opens  the 
question,  whether  the  characteristic  length  is  redundant  altogether  since  distributed 
stress  concentrators  at  aggregate  vertices  trigger  crack  initiation  in  a  fairly  distributed 
fashion? 

We  recall  that  associated  flow  results  in  highly  dilatant  behavior  and  in  excessive 
confinement  when  plane  strain  conditions  are  adopted.  In  return,  this  passive  con¬ 
finement  leads  to  gross  overprediction  of  mortar  strength  when  high  dilatancy  angles 
are  used.  Thus  plane  stress  conditions  are  adopted  for  the  inelastic  response  analysis 
in  which  the  internal  friction  angle  of  $  =  30®  assures  closure  of  the  failure  envelope 
in  biaxial  compression.  The  analysis  is  carried  out  with  the  nonlinear  finite  element 
program  Microfem  [32]  in  which  the  invariant  response  formulation  [33]  has  been  im¬ 
plemented  for  the  general  description  of  frictional  materials.  Fig. 11  summarizes  the 
properties  of  the  Drucker- Prager  model  and  the  softening  properties  in  uniaxial  ten- 
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sion.  Although  the  low  softening  modulus  results  in  far  too  ductile  behavior  in  direct 
tension,  it  is  requited  to  assure  stable  post-peak  control  and  to  avoid  snap-through 
phenomena  particularly  in  biaxial  compression. 

The  material  behavior  of  the  mortar-aggregate  interface  is  depicted  in  Fig.  12  for 
tensile  failure  which  is  mobilized  when  the  current  value  of  degrading  bond  strength  is 
exceeded.  This  approach  monitors  constant  values  of  the  fracture  energy  release  rate  in 
tarms  of  the  crack  opening  displacement,  whereby  the  shear  stiffness  diminishes  to  zero 
when  the  residual  tensile  strength  has  been  exhausted.  In  sum,  the  present  interface 
model  allows  to  investigate  two  extreme  cases,  (i)  full  bond  without  degradation  of 
the  interface  behavior,  and  (ii)  purely  tensile  debonding  without  shear  transfer  beyond 
critical  crack-opening. 

4  Computational  Results 

The  investigation  of  the  matrix-aggregate  composite  is  comprised  of  four  parts:  The 
first  study  is  focussed  on  the  elastic  properties  of  the  composite  for  various  ratios  of 
elastic  phase  moduli.  The  next  two  investigations  examine  the  tensile  and  compressive 
failure  mechanisms  in  rectangular  composite  bodies  with  different  aspect  ratios.  The 
fourth  study  considers  compressive  and  tensile  load-histories  during  which  a  composite 
specimen  is  preloaded  in  one  direction  and  subsequently  unloaded  and  reloaded  in  the 
orthogonal  direction  up  to  failure.  In  most  cases  monolithic  contact  i.e.  full  bond  is 
assumed  between  aggregate  and  mortar.  In  few  cases  tensile  debonding  is  permitted 
in  order  to  study  the  influence  of  degrading  interface  properties  surrounding  each 
aggregate  particle  during  failure  progression. 

4.1  Effective  Elastic  Properties  of  Particle  Composite 

In  the  first  study  we  examine  the  effect  of  different  elastic  mortar  and  aggregate  prop¬ 
erties.  This  topic  is  normally  pursued  within  the  framework  of  composite  mechanics 
with  mixture  rules  which  define  the  effective  stiffness  properties  of  the  overall  com¬ 
posite.  The  results  of  the  two-phase  particle  composite  are  shown  in  Fig.  13  for  the 
random  distribution  of  polygonal  inclusions.  The  four  diagrams  depict  the  change  of 
composite  stiffness  when  the  ratio  between  aggregate/matrix  moduli  varies  between 
one  and  ten  for  the  aggregate  concentration  of  50%. 

Clearly,  the  random  arrangement  of  inclusions  leads  to  composite  stiffness  proper¬ 
ties  which  lie  between  the  two  limiting  cases  of  the  series  and  parallel  arrangement  of 
the  two  constituents.  The  results  of  the  numerical  simulation  agree  surprisingly  well 
with  the  analytic  estimates  by  Hashin  [34]  for  non-contiguous  spherical  inclusions  and 
low  level  particle  concentration.  Altogether  the  graphs  of  Ftp.  13  indicate  that  the 
series  model  provides  a  lower  bound  of  the  overall  stiffness  which  exceeds  the  mortar 
stiffness  merely  by  70%.  On  the  other  hand,  the  parallel  model  leads  to  a  lir.jar  in¬ 
crease  of  the  composite  stiffness  which  grows  proportionally  to  the  aggregate  stiffness 
and  the  phase  concentration. 


The  angular  particle  composite  exhibits  elastic  stiffness  properties  which  are  asymp¬ 
totically  controlled  by  the  elastic  properties  of  the  contiguous  matrix  surrounding  each 
particle.  In  fact,  the  normalized  values  of  the  effective  composite  stiffness  inreases  in 
our  case  not  more  than  three  times  the  value  of  the  matrix  stiffness.  The  same  obser¬ 
vation  holds  for  the  composite  strength  in  tension  which  is  essentially  controlled  by 
that  of  the  mortar  matrix  with  //  =  0.82  Ksi,  see  /’fj.14.  This  figure  also  illustrates 
the  isotropy  of  stiffness  and  strength,  when  the  composite  body  is  subjected  to  two 
separate  tensile  load  histories  in  x-  and  y-directions.  In  contrast,  Fig. 15  depicts  the 
results  of  two  analogous  experiments  in  uniaxial  compression.  Again,  the  stiffness  of 
the  composite  body  exhibits  full  isotropy  as  expected,  however  the  uniaxial  strength 
values  differ  considerably  in  x-  and  y-directions.  Note  that  both  composite  strength 
values  exceed  the  nominal  strength  f\  =  2.0  Ksi  of  the  mortar  matrix,  in  spite  of  the 
homogeneous  boundary  conditions  adopted  for  the  numerical  experiments.  In  con¬ 
trast  to  the  tensile  strength  of  particle  composites,  the  heterogeneous  microstructure 
increases  the  compressive  strength  of  the  composite  body  beyond  that  of  the  mortar 
matrix,  because  of  the  locking  block  action  of  the  aggregate  particles  embedded  in  the 
frictional  mortar  matrix. 

Three  specimen  geometries  are  examined  in  order  to  address  the  issue  of  surface 
versus  volume-dominated  failure  processes  in  direct  tension  and  uniaxial  compression. 
The  results  of  the  square  specimen  with  the  aspect  ratio  of  h/l  =  1  are  denoted  by  the 
symbol  (4*4),  while  those  of  the  rectangular  specimens  with  the  aspect  ratios  h/l  =  2 
and  h/l  =  0.5  are  designated  by  the  symbols  (8  *  4)  and  (2*4).  The  stiffness  ratio 
Ea/ Em  —  5  is  used  throughout  the  investigation  for  the  matrix-aggregate  composite. 
The  gradually  increasing  load  is  applied  under  mixed  control  with  full  lateral  restraint 
at  the  top  and  bottom  faces  of  the  specimen  simulating  the  effect  of  rigid  loading 
platens.  Thus,  both  the  micro-structure  as  well  as  the  lateral  restraint  influence  the 
strength  and  the  localization  process  at  failure. 

4.2  Tensile  Failure  Simulations 

Fig.16  illustrates  the  overall  response  behavior  of  the  square  specimen  subjected  to  di¬ 
rect  tension.  The  results  are  plotted  in  terms  of  nominal  values  of  stresses  and  strains 
which  are  extracted  from  the  reactions  and  the  relative  displacements  between  the 
boundaries  of  the  finite  element  idealization. The  response  predictions  compare  well 
with  experimental  results  in  direct  tension,  particularly  when  the  localization  at  peak 
and  the  progressive  separation  in  the  post-peak  regime  are  considered.  The  active 
zones  of  plastic  dissipation  are  inserted  in  Fig.  16  in  order  to  visualize  the  progres¬ 
sive  failure  process  in  tension.  The  three  plots  exhibit  different  failure  stages  in  the 
pre-  and  post-peak  regimes.  The  distributed  cracking  mode  in  the  pre-peak  regime 
is  followed  by  strong  redistribution  during  the  formation  of  a  single  transverse  crack 
which  separates  the  composite  specimen  from  the  loading  platen  near  the  bottom  face. 

The  post-peak  results  of  the  three  specimens  with  different  heights  are  summarized 
in  Fig.ll.  The  tall  specimen  exhibits  only  slightly  more  brittle  behavior  during  sep- 
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aration  than  t.ie  compact  specimen.  In  fact,  the  three  post-peak  slopes  do  not  differ 
enough  to  contradict  the  concepts  of  fracture  mechanics  for  interpreting  the  separa¬ 
tion  process  in  tension.  Recall,  that  the  area  under  the  nominal  tensile  stress  versus 
c*-ack  opening  displacement  curve  is  a  measure  of  the  fracture  energy  release  rates  for 
the  three  different  specimen  heights.  The  agreement  confirms  previous  experimental 
observations  on  tensile  failure  of  different  height  specimens  [35],  where  the  fracture 
energy  release  rate  remained  fairly  constant  during  formation  of  the  transverse  crack. 
In  sum,  tensile  softening  is  a  surface-  rather  than  a  volume-controlled  fracture  phe¬ 
nomenon. 

The  question  of  the  characteristic  length  is  addressed  in  .Ft0.I8.  The  plot  depicts 
the  locations  of  interface  elements  in  the  square  specimen  which  are  subjected  to  tensile 
failure  at  an  overall  load  level  corresponding  to  75%  of  the  tensile  peak  load.  Recall 
that  the  polygonal  particles  are  surrounded  by  interface  boundary  layers  with  50% 
higher  debonding  strength  than  the  contiguous  mortar.  This  premature  debonding  in 
the  interface  elements  is  mobilized  in  spite  of  their  high  debonding  strength,  because 
the  interface  stresses  are  dominated  by  the  high  stress  level  in  the  stiff  aggregate  par¬ 
ticles.  In  spite  of  these  stress  concentrators  the  tensile  cracking  path  is  ultimately 
forced  through  the  mortar  matrix  because  of  the  non-contiguous  distribution  of  aggre¬ 
gate  particles  and  the  disconnected  network  of  interface  boundaries.  The  plot  clearly 
underlines  the  distributed  nature  at  crack  initiation  which  is  governed  by  the  loca¬ 
tion  of  aggregate  particles.  However,  the  final  separation  process  develops  in  the  form 
of  a  single  discontinuity  across  the  mortar  matrix  with  contiguous  zones  of  elastic 
unloading.  Thus,  we  may  conclude  that  the  final  separation  process  is  a  structural 
phenomenon  which  affects  the  entire  specimen  due  to  localization  into  active  zones  of 
high  energy  dissipation  and  large  contiguous  blocks  of  elastic  unloading. 

4.3  Compressive  Failure  Simulations 

For  tensile  strength  degradation  it  is  widely  accepted  that  tensile  cracking  is  governed 
by  fracture  mode  I.  The  question,  whether  fracture  mechanics  also  governs  the  fail¬ 
ure  process  in  compression,  is  debated  vigorously  with  respect  to  the  issue,  whether 
there  exists  a  critical  fracture  energy  release  for  mode  II-type  fracture  for  brittle  solids 
independently  of  mode  I.  For  this  reason  the  numerical  simulations  of  the  tension  ex¬ 
periments  above  are  now  repeated  for  the  case  of  uniaxial  compression. 

The  compressive  response  behavior  is  shown  in  Fig.  19  for  the  square  composite 
specimen.  The  inserts  in  the  figure  provide  some  insight  into  the  failure  patterns 
which  develop  in  the  pre-  and  post-peak  regimes.  The  elementary  Drueker-Prager 
model  of  the  mortar  matrix  does  result  in  shear-type  failure  modes  which  agree  sur¬ 
prisingly  well  with  the  experimental  observations  reported  by  von  Mier  [7].  Although 
the  multiple  shearing  planes  are  partially  induced  by  the  boundary  conditions,  the  final 
conical  failure  mode  emerges  along  weak  shear  planes  in  the  contiguous  mortar  matrix. 
Opposite  to  the  response  behavior  in  direct  tension  the  matrix-dominated  shear  failure 
leads  to  an  increase  of  the  uniaxial  compressive  strength  of  the  composite  beyond  the 
mortar  strength.  The  numerical  simulation  illustrates  the  frictional  strengthening  due 
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to  composite  action.  In  that  case,  both  strength  and  softening  behavior  are  clearly 
influenced  by  the  microstructure  of  the  mortar-aggregate  composite. 

The  issue  of  surface-  versus  volume-dominated  failure  mechanisms  in  compression 
is  addressed  in  Fig.20.  The  plots  depict  the  post-peak  response  behavior  of  three 
composite  specimens  with  different  heights.  The  different  softening  slopes  indicate 
that  the  fracture  energy  release  rate  varies  strongly  with  specimen  height,  quite  in 
contrast  to  the  experimental  observations  reported  by  von  Mier  [7].  On  the  other 
hand  it  is  intriguing  that  normalizing  the  axial  displacements  by  the  three  specimen 
heights  leads  to  stiain-softening  response  diagrams  which  agree  fairly  closely.  In  other 
terms,  the  numerical  experiments  in  compression  support  the  notion  of  volume-  rather 
than  surface-dominated  failure  processes. 

4.4  Compressive  and  Tensile  Failure  Simulations  with  Or¬ 
thogonal  Preloading 

In  the  final  study  we  examine  the  effect  of  preloading  on  the  stiffness  and  strength 
of  square  composite  specimens  in  which  the  direction  of  loading  is  rotated  by  90°  at 
different  states  of  damage.  To  this  end  a  number  of  virgin  composite  specimens  were 
first  loaded  uniaxially  in  the  y-direction  up  to  a  prescribed  level  of  preloading  and 
subsequently  unloaded  and  reloaded  in  the  x-direction. 

The  compressive  results  are  shown  in  Fig.2\  depicting  the  strong  degradation  of 
compressive  strength  because  of  preloading  in  the  orthogonal  direction.  The  analogous 
simulation  in  tension  shown  in  Fig. 22  exhibits  far  less  strength  degradation  inspite  of 
tensile  preloading  beyond  peak.  The  same  preloading  effects  in  tension  and  compres¬ 
sion  were  reported  by  von  Mier  [7]  who  discerned  in  his  concrete  experiments  between 
‘cylindrical’  and  ‘planar’  modes  of  failure  which  are  accompanied  by  totally  different 
changes  of  the  internal  microstructure.  In  compression,  the  inclined  mode  of  failure 
takes  place  simultaneously  in  two  orthogonal  directions.  Therefore,  reloading  in  the 
orthogonal  direction  is  strongly  influenced  by  the  level  of  preloading. 

Fig.22  illustrates  the  different  stages  of  inclined  failure  during  compressive  preload¬ 
ing  in  the  y-direction  and  subsequent  compressive  loading  in  the  x-direction.  In  con¬ 
trast,  Fig. 24  depicts  the  localized  failure  modes  in  tension  which  are  virtually  orthog¬ 
onal  to  the  two  tensile  load  directions  with  little  interaction. 

In  sum,  the  numerical  load  history  experiments  clearly  exhibit  the  directionality  of 
progressive  failure  and  the  related  differences  between  tensile  and  compressive  strength 
degradation  of  specimens  preloaded  previously  in  orthogonal  direction. 

5  Concluding  Remarks 

The  composite  simulation  provides  new  insight  into  the  failure  mechanisms  which  take 
place  in  the  heterogeneous  matrix-aggregate  composite.  The  computational  results 
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indicate  that  the  initial  failure  mechanisms  are  very  distributed  both  in  tension  and 
compression.  Only  the  final  stage  of  the  tensile  separation  process  is  truly  localized. 
In  that  case  fracture  mechanics  prevails  in  form  of  a  surface-dominated  process  with  a 
characteristic  length  which  is  of  the  same  order  of  magnitude  as  the  specimen  geometry. 
This  localization  mode  is  accompanied  by  dramatic  strain  redistribution  which  requires 
dynamic  enrichment  or  realignement  of  the  finite  element  mesh  and  endowment  of  the 
constitutive  operator  by  a  characteristic  length  to  delimit  the  zones  of  active  failure. 
On  the  other  hand,  the  failure  mechanisms  in  compression  exhibit  fracture  properties 
which  differ  substantially  from  those  in  tension. 
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Fig.  1  Direct  Tension  Experiments  by  B.  Hurlbut[6j 
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Fig.  2  Triaxiai  Compression  Experiments  by  S.  Smith[8] 


Fig.  3  Failure  Modes  in  Triaxiai  Compression  and  Extension 
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Fig.  4  Incremental  Plastic  Strains  Superposed  on  Triaxial  Stress  Path 
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Fig.  5  Wedge  Experiment:  Stress  Path  and  Incremental  Plastic  Strain  Vectors 
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Fig.  6  Incremental  Plastic  Work  at  Different  Stages  of  Wedge  Experiment 
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Triangulation  of  Mortar  Phase 
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Fig.  11  Decohesive  Drucker-Prager  Description  of  Mortar  Phase 


Fig.  12  Interface  Element  for  Tensile  Debonding  of  Mortar-Aggregate  Contact 
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Fig.  18  Tensile  Failure  Pattern  of  Interface  Debonding  at  75%  Peak  Strength 
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Fig.  14  Isotropy  of  Stiffness  and  Strength  in  Uniaxial  Tension 
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Fig.  15  Isotropy  of  Stiffness  and  Strength  in  Uniaxial  Compression 
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Fig.  16.  Progressive  Failure  Behavior  in  Direct  Tension  Experiment 


Fig.  17  Post-Peak  Behavior  of  Different  Height  Tension  Specimens 
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Fig.  20  Post-Peak  Behavior  of  Different  Height  Compression  Specimens 
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Fig.  23  Summary  of  Load  History  Experiments  in  Compression 
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Fig.  24  Summary  of  Load  History  Experiments  in  Tension 


